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Abstract 

The anomalous dimensions of four-quark operators (qiqj)v—A(QkQl)v-A are calculated 
in the large Nf limit. As expected, the result is a convergent series without renormalon 
ambiguities. Using the approximation of "Naive Nonabelianization" , an additional all- 
order contribution to the anomalous dimension matrix is obtained which is somewhat 
larger than the exact NLO correction itself. Possible phenomenological applications in 
nonlcptonic B decays and B — B mixing are briefly considered. 
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1 Introduction 



The behaviour of the perturbative series at large orders constitutes one of the major unsolved 
problems of quantum field theory. Since it is in practice impossible to calculate a given 
physical quantity exactly to all orders of perturbation theory, the following approach for a 
better understanding of the a priori unkown large-order behaviour suggests itself: namely, 
instead of considering all possible diagrams, one tries to identify certain subclasses of diagrams 
which become dominant at sufficiently high orders and can be calculated exactly to all orders. 

Once such a subclass of diagrams is found, it will be useful in three aspects: (i) for analyz- 
ing fundamental questions related to the summability and/or uniqeness of the perturbative 
series, as well as the connection of these problems to non-perturbative issues, (ii) for estimating 
the numerical effects of higher order contributions, which are important in many phenomeno- 
logical applications, (iii) for having at one's disposal an independent check of forthcoming 
exact higher-order calculations. The latter point should especially be emphasized, since such 
calculations will most certainly be performed in a completely computerized way. Therefore, 
analytical expressions for a well-defined part of the full result should be regarded as a helpful 
and desirable tool for the testing of these probably very extensive computer programs. 

In fact, until now only one of these subclasses of Feynman graphs has been found : di- 
agrams with an insertion of an arbitrary number of massless fermion loops ("bubbles") in 
a gluon line, the so-called renormalon chains. If a physical quantity R has the perturbative 
expansion 

oo 

R=Y1 r n" n , (1) 
n=0 

where a is the coupling constant of the theory, and if we perform for each coefficient r n an 
additional expansion in Nf, the number of massless fermions in the theory, 

n-1 

r n = £ r ^ N T^ n>l, (2) 

m=0 

then the coefficients rl n ^ will obviously be determined by diagrams where n — 1 fermion 
bubbles have been inserted into a gluon line of some leading-order graph. 

Of course it is doubtful whether these coefficients ri™ ^ correctly describe the asymptotic 
behaviour of the perturbative series or not. Although this will be the case in the limit 
Nf — > oo, at first sight it looks like a gross mistake to conclude that the real world (with 
Nf = 5 in the case of B physics or e + e _ annihilation) resembles at least in some aspects this 
fictitious limit. Yet there are some arguments for such a point of view. Above all they rely 
on the observation ffl,H that the amplitude of an individual fermion-bubble diagram grows 
in general like n\ at high orders. If one remembers that usually not the amplitude of an 
individual graph but only the number of all possible graphs grows with n\, and if one further 
assumes that all diagrams yield approximately the same contribution with no preferred sign, 
then one can argue that the fermion-bubble diagrams should become dominant at sufficiently 
high orders. For recent reviews and further references see ||, |J. 
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In practical applications, higher order effects play a significant role only in QCD with its 
relative large coupling constant. However, QCD is exceptionally bad described by the limit 
Nf — > oo, since in this limit asymptotic freedom, one of the most characteristic features of the 
strong interaction, is violated (/3q ~ Nf — 33/2, the first coefficient of the QCD beta function, 
becomes positive for Nf > 16). As a generalization of the QED case, where (3q ~ Nf, it was 
therefore proposed ^-0] that in QCD one should expand the perturbative coefficients r n in 
powers of Nf — 33/2, 



n-1 

\ " Wri 



r n = 2^rW(N f -33/2) m , n > 1, (3) 

and that physical quantities are at sufficiently high orders well described by neglecting all but 
the coefficients rjf ^ = ri™ ^ in this expansion. This approximation is usually referred to as 
"Naive Nonabelianization" (NNA). 

In the last two years the NNA method was applied to several observables where a compari- 
son with existing Next-To-Leading-Order (NLO) or Next-To-Next-To-Leading-Order (NNLO) 
calculations was possible. Specifically, analyses were performed of the R ratio and related 
quantities J5J , of the hadronic decay width R T of the r lepton Q , of semileptonic B decays § , 
and of structure functions in deep inelastic scattering |i~(i|-13|. In all of these, a rather good 



agreement between NNA and the exact calculation was claimed at third or even at second 
order. However, it should be stressed that a systematic reason for this agreement is still not 
known. 

It is interesting to extend this program to some other quantities were exact NLO results 
are available. To this end, in this paper the anomalous dimensions of the weak current-current 
four-quark operators is calculated in the large Nf limit. Such operators occur in the effective 
Hamiltionian for e. g. weak hadronic decays and particle-antiparticle-mixing involving mesons. 
There NLO anomalous dimension matrix was calculated in [14], Dj, see also [16] for a recent 



review on the subject of effective hamiltonians for weak decays. Clearly, those anomalous 
dimensions are no direct physical observables; a complete calculation of the effective hamilto- 
nian additionally includes the matching of the effective onto the full theory at the electroweak 
scale. This procedure involves the calculation of finite parts of Feynman diagrams. At a given 
order in a s , the matching is therefore considerably more difficult than the computation of 
renormalization group functions which is performed in what follows; nevertheless a posteriori 
it turns out that it is presumably even more important. This point will be discussed in the 
last but one section. 

The organization of this paper is as follows: In section 2, the scene is set with the basic 
definitions of the relevant operators and their anomalous dimensions; it is then explained how 
to obtain these anomalous dimensions from the calculation of truncated Green functions. 
In section 3, this calculation is performed and the necessary renormalization of the gluon 
field in the large Nf limit is discussed in detail. From this an analytic large Nf formula 
for the anomalous dimensions, summed to all orders in a s , can be obtained. This is done 
in section 4. The result is discussed in section 5, where the NNA approximation is used to 
estimate the perturbative contributions beyond next-to-leading order, an estimate that will 
be relevant for some phenomenological applications. It is also discussed how the result for 
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the anomalous dimensions should be seen in the context of the full large Nf calculation of 
the effective Hamiltonian. The paper closes with a short summary in section 6. Finally, an 
appendix is devoted to the so-called evanescent operators, objects that occur in the course of 
the calculation in section 3. In particular, it is proved that these operators do not contribute 
to anomalous dimensions in the large Nf limit. 



2 Definitions 

Current-current operators arise due to tree-level W exchange. There Dirac structure is given 
by 

7 = %7 M (1 - 7s)?i %7 M (1 - 7s)«, (4) 

where k, I denote flavour indices. Due to exchange of gluons between the quark legs, two 
different color structures can arise, in symbolic form written as 

1 = <W <5 7 <5, 1 = 5 a s 5pj, (5) 

where a, /3, 7, 5 denote the color indices of the first, second, third and fourth quark field q in 
Eq. (Q). In this compact notation, the two operators we are dealing with are given by 

O x = 7 1, 

2 = 7 1. (6) 
For further purposes we also introduce the linear combinations 

± = l -{0 2 ±0 1 ). (7) 

In this basis, the anomalous dimension matrix will be diagonal. If we now calculate some 
truncated Green function (Oj) with an insertion of one of these operators, the result will 
in general be divergent even after renormalization of fields and QCD coupling. In order to 
achieve a finite result, we make use of our freedom to redefine the operator basis0 according 
to 

Oi = ZijOf. (8) 
If we additionally renormalize the four external fields, 

Q = zl /2 q R , (9) 
then we can write the truncated Green function as 

(O i ) = Zj 2 Z ij (Of), (10) 



1 - ) In general, one has to include additional "evanescent" operators in the operator basis, since they unavoid- 
ably arise while calculating loop diagrams like those of Fig. 1. Accordingly, these evanescent operators also 
participate in the multiplicative renormalization of Eq. (B). However, in the large Nf limit they can safely be 
neglected, as shown in the appendix. 
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where the (O^'s are now finite. The subtractions will be performed in the MS scheme 
throughout, so (working in D = 4 — 2e dimensions) the renormalization constants Zij will 
have the form 

= ! + £-£- (11) 

k=l 

Zij is a function of the renormalized coupling and therefore depends on the renormalization 
scale \x. The anomalous dimension matrix of the operators Oi is defined as 

l l3 = {Z-% klJ ^-. (12) 

Introducing the quantity 

Zij = Z^Zij, (13) 



which can, according to Eq. (10), be read off directly from the divergent pieces of the 
truncated Green functions (Oi), we may write Eq. (12) as 

7y = lij + 2 7</>, (14) 

with 

~ i dZfcj i dZyj] 

ji j = (z-) ik »—, li} = z;»—. (15) 

Using the equivalents of Eq. ( JTl] ) for Zij and Z$, the chain rule d/dfj, = (dg/dfi) d/dg (where g 
is the renormalized QCD coupling) and the renormalization group equation for g one obtains 



% = - 2g2 ^-' ^ = - 25 ^' (16) 

i.e. the anomalous dimensions are directly related to the 1/e pole part of the corresponding 
renormalization constant. 



3 Calculation of Green functions with bubble insertions 

We now proceed to calculate the truncated Green functions (0\) and (O2) with an insertion 
of n fermion bubbles into the gluon line, i.e. the diagrams shown in Fig. 1. In order to 
determine the quark wave function anomalous dimension 7^, we will additionally calculate 
the truncated two point function depicted in Fig. 2. To this end, we need as a first step 
an expression for the gluon propagator with an insertion of n > fermion bubbles (but no 
counterterms) . A straight forward calculation in D = 4 — 2e dimension yields 

^(k) = ^Uv-^)\-m 2 )] n - (17) 
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Figure 1: The six Feynman graphs relevant for the anomalous dimensions of the operators 0\ and 
02- The dashed line denotes the gluon propagator with an insertion of n fermion bubbles and the 
appropriate counterterms. 




Figure 2: The Feynman graph relevant for the anomalous dimension of the quark wave function. 
The dashed line again denotes the gluon propagator with an insertion of n fermion bubbles and the 
appropriate counterterms. 



Here, the function U(k 2 ) stems from the one loop calculation of the fermion bubble, 



U{k 2 



a.- 



N 



7T 



/ 



-k 2 



r(e)B(2-e,2-e) 



(18) 



and k is the momentum flowing through the propagator. It should be stressed that Eq. (|lj 
is valid in any gauge for n > 1. Only for n = it is the specific result of the Landau gauge. 
So the fermion bubble contribution is for itself a gauge independent quantity, as of course 
must be the case. 

Now one can insert the effective propagator (|i~7|) into the diagrams of Fig. 1 and per- 
form the loop integrations via introduction of Feynman parameters as usual. Adding the 
contributions from all six diagrams, we find 

\ bare _ 

' ■ F(e, (n + l)e) M c 



{Ox 



(O, 



\ bare 



47T 



a s N } 
6tt 



■1) 



ip^tree 



(n + l)e n+1 



with 



F(e,z) 



47T/X 



r(l + z)r(l-z) 6er(e)B(2 - e,2 - e) 



<0s 



z/e-i 6-4e 



\tree 



, (19) 



r(3 



(20) 
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1 + e 

and Cf = (N 2 — 1)/(2N), N = 3. (Oi) tree denote the trivial tree-level matrix elements of the 
operators Cj, and the auxiliary function F(e, z) has been defined in such a way such F(0,0), 
F(e, 0) and F(0, z) are finite. From a calculational point of view, Eqs. (|l9| ) - (f2~l|) constitute 
already the main result of this paper, so several remarks are in order: 

(i) The index "bare" indicates that until now no renormalization has been performed. In 
the next step we will renormalize the gluon field so that the transition (Oi) bare — > (0{) is 
achieved. The final operator renormalization (Oi) — ► (Cf 2 ), cf. Eq. (||), is the subject of the 
following section. 

(ii) In calculating Eq. (|T^) we set all external momenta to zero which is allowed since the 
renormalization constants we are interested in do not depend on the external states. Infrared 
divergent integrals were regulated by introducing a finite gluon mass A that appears as the 
only scale (apart from the renormalization scale fj) in Eq. (pC|). 

(iii) In order to obtain Eq. (|l9|) we had to project complicated Dirac structures, generically 
of the form r ® T, onto the simpler V — A (g) V — A structure of the operators Oi, i. e. we used 
identities like 

yyy (1 _ 75 ) ® llxlvlx (l - 75 ) = 4 (4 - e - e 2 ) -f{l - 75) ® 7^(1 " 7s)- (22) 

A general method for obtaining such relations is explained e. g. in |18j . However, such a 
projection is a very peculiar procedure since it has to be done in D = 4 — 2e dimensions 
where, strictly speaking, no complete finite basis {7^ ®> 7^} of Dirac structures on which 



to project can be given. Eq. (22) should therefore be understood as the projection on a 
certain subset of such a basis, namely, the one that forms the complete basis in D = 4 
dimensions. Only in this sense the coefficient 4 (4 — e — e 2 ) in Eq. (^) is a well defined 
quantity. Rigorously, the discrepancy between the l.h.s. and the r.h.s. of Eq. (|22| ) defines 
an evanescent operator (vanishing in the limit e — * 0) which has to be treated correctly in 
higher order calculations ||l8| . However, it is shown in the appendix that such operators are 
irrelevant for the purpose of this paper. 

In order to renormalize the gluon field, the quite ingenious method described by Beneke 
and Braun in the appendix of [19] is used. One first needs the counterterm appropriate to 



cancel the divergence of one fermion bubble. From Eq. (jl8|) it can be read off to be 
in the MS scheme. Since there are n!/(fc! (n — k)\) possibilities to replace k out of n fermion 
bubbles with this counterterm, and since all these possibilities have to be added, one can infer 
from Eq. ( |i~9l) the following expression for the (partially) renormalized Green functions: 

1 






'a a N f \ 


4vr 1 


v 6vr J 


/ n 


) ("l) fc 


[ k 


/ fc+1 



F(e, (* + !)€) 




{Oi) tree 
(0 2 ) tree 



(23) 



We called this the partially renormalized Green functions, because operator and wave function 
renormalization are still to be performed. 
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It is gratifying that the finite sum in Eq. (|23|) can be solved. Expanding F(e, z) in powers 



of z, 



j=0 



we have 

n 



E 

k=0 



n 
k 



(~l) fc 



F(e, (fc + l)e) 



n+1 
j=0 



E 

.fc=0 



n 



(_l)fe( fc + 



(24) 



(25) 



The sum over j stops at j = n + 1, since the remaining terms do not contribute in the limit 

e — ► (note that Eq. ( p5| ) is still to be multiplied with an overall factor l/e n+1 ). Using 

appropriate formulas from [2C,21], the sum over k in Eq. ( f25| ) can easily be performed. It is 

different from zero only in the cases j = and j = n + 1, 

\& i 

/ (e)_ _ + / n+1 ( e ) e »+i(-l)»n! . (26) 



E 

fc=0 



A; 



-1) 



F(e, (fc + l)e) 



n+l 



The term proportional to e n+1 can be discarded since it yields only a finite contribution in 
the limit e — > which is irrelevant for the anomalous dimensions. With 

3- 2e 



Me) 



3r(l + e)B(2 - e, 2 - e)T(3 - e) 
from Eq. (^0|), one finally obtains 



(27) 



<Oi) 



47T 



3- 2e 



1 



(n + l)e n+1 



3r(l + e)B(2 - e, 2 - e)T(3 - e) 



+ finite 



M r 



(01) tree 

(0 2 ) tree 



(28) 



From this formula the renormalization constants Zij can immediately be read off. 

A completly analogous calculation can be done for the fermion self energy with an insertion 
of n fermion bubbles, see Fig. 2. This calculation yields 



i£ (n) (p) 



%i> U + 



47T 
1 



6tt 



e(3 - 2e)r(4 - 2e) 



+ finite 



(29) 



.(n + l)e™ +1 6r(l + e)r(2 - e) 2 r(3 - e) 

from which the wave function renormalization constant can be read off. Eq. ( p9[ ) is again 
gauge independent for n > 1 and specific to Landau gauge in the case n = 0. 



4 Extraction of anomalous dimensions 

Summarizing the results of the previous section, renormalization constants in the large Nj 
limit may in general be written as 

^a s (asN f \ k f (e) 
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at any fixed order n + 1 in perturbation theory. Here, /o(e) is some function of e which may 
also contain color factors. However, we prefer to have at our disposal an all-order result for 
the corresponding anomalous dimension. This can be achievied by expanding /o(e) in powers 
of e, 

00 

/ ( e ) = £/; ei - (31) 

Inserting this into Eq. (|30|), we obtain for the 1/e part of the renormalization constant, in the 
notation of Eq. (11), 

z(1) = _^M^i (32) 
^ 4vr V 6vr J k + 1 ' y ' 

which can formally be summed up to all orders by setting n to infinity, 

j =0 Att{ 6tt ) k + l 
From this the corresponding anomalous dimension is obtained according to Eq. ( |l~6| ) as 

dZ (l) 

n all— orders 

1 all— orders — ^0L S 



7 (i) _ ^ a s f a s N f \ k f k 

^all-orders ~ L ^ ft, Ixl" 



da s 

k=0 



2ir V 67r 



^Aliri- (34) 



where in the last step Eq. fl3l| ) was employed. 

Thus we get from Eq. (^) for the anomalous dimension of the quark field 

^ = ^C F fJ^l), (35) 



2-7T \ 6tt 

with 

„ , , 1 . . r(4-2x) 

= --x(3 - 2 -) r(1 + x)r(2 _ x)2r(3 _ x) ' (36) 



in agreement with a result that Gracey |22j derived some years ago using a completely different 
method. Having found this, we finally obtain from Eqs. (|14|) and ( ]28| ) the following formula 
for the anomalous dimensions of the operators 0\ and O2 in the large Nf limit: 

w -* + ^-£, (=^)(-f _/v , 37 , 



with 



1 IY4 — 2x) 

fo W = 5 (i + »)(»- a») r(1 + ^ _ I)2 ' r(3 - ,) • < 38 > 
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This can alternatively be written in the diagonal operator basis (pi), where one has 

7± = ± ^ M^rJ ~- (39) 

Reexpanding in powers of a s , our result (|39| ) agrees with the atgNf part of the exact NLO 
result as given in Eq. (5.2) of [jlSf] - One should mention that although our findings were derived 



in the MS scheme, they are of course also valid in the MS scheme (with the replacement 

a MS ^ a MS 

understood). This is because the MS scheme is related to the latter simply 
by a multiplicative scale redefinition /i — ► (e 7B /47r) 1//2 /i, so that the functional dependencies 
of renormalization group functions on a s do not change. On the other hand, it should be 
emphasized that the result for the anomalous dimensions j± in general as well as in the large 
Nf limit depends on the regularization scheme used for 75 in D dimensions: all our results 
were obtained in the "Naive Dimensional Regularization"' (NDR) scheme with a naively 
anticommuting 75. 



5 Discussion of results 

As explained in the introduction, the large Nf limit of QCD unfortunately corresponds by 
no means to the physical world. Nevertheless it is widely conjectured that one can still make 
some contact to reality, if one uses the results of the large Nf limit as a starting point for 
the so-called NNA approximation, i. e. if one replaces Nf by Nf — 33/2 in these results. This 
procedure, motivated by the analogy to QED where Nf is essentially the first coefficient of the 
QED beta function, may also be called the large (3q approximation. With (3$ = 2Nf/3 — 11, 
it amounts to the replacement 



a s Nf a s p 



(40) 



6ir 4tt 

in Eq. (j39|). To the knowledge of the author, the only justification for this procedure is a 
certain phenomenological success; no further comments on its validity will therefore be found 
in this paper. Instead, the substitution ( [40|) will be used as an admittedly crude and naive 
prescription for estimating higher-order effects. 

Before doing so it is interesting to investigate how well the NNA approximation works 
already at NLO, where a comparison with the exact result is possible. Expanding 

(0) ot s (i) ( a s x -' 



7± = 7r£+7 w £ ) +..., (41) 



the exact result [15| in the NDR scheme reads 



+ 7« = -^-f. (42) 

In the NNA approach, the constant terms are neglected. With Nf = 5 or /3q = —23/3, 
this is obviously a very good (accurate to 7%) approximation for 7+ , but also a very bad 
(wrong sign) approximation for 7_ . So one cannot claim from this comparison that the NNA 
approximation generally works well at low orders^. 

2 'See, however, the remark at the end of this section. 
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Despite this somewhat discouraging observation, let us now ask for the contributions of 
all orders beyond NLO in the NNA approach. Evaluating the anomalous dimensions at a 
scale n ~ 5GeV or a s (fi) ~ 0.21, we have 

7+ = +0.0670 (LO) - 0.0014 (exact NLO) - 0.0021, 

7 _ = -0.1350 (LO) - 0.0053 (exact NLO) + 0.0043. (43) 

Here, the first two numbers are exact result and the third number indicates the summed 
contribution from all orders beyond NLO in the NNA approximation (i.e. the terms o^0q, 



a^0Q, . . .), obtained from Eq. (|39|). We notice that this contribution is of the same order of 
magnitude as the exact NLO contribution itself. 

One can go one step further and examine to what extent the Wilson coefficients (//) of the 
operators Oi are affected by these additional fermion-bubble contributions to the anomalous 
dimensions. We will investigate here two cases that are of some phenomenological interest: 
nonleptonic B decays and B° — B° mixing. 

(i) Nonleptonic B decays 

Only decays of the b quark into three different flavours, e. g. b — > cud, are considered. For 
these decay modes, no penguin operators can occur and the complete effective Hamiltonian 
is given by 

Gf 
72 

with the four-quark operators Oi as defined in Eq. (ph. The renormalization group analysis 
is most easily be done in the diagonal basis 0± with the corresponding Wilson coefficients 
z± = C% +_ C\. These Wilson coefficients obey the renormalization group equation 

/^TJ^fcOO = 7±z±0«) (45) 



H eff = -4 Vckm [Ci(/i)0i + C 2 (fJ,)0 2 ], (44) 



which can be solved numerically using the anomalous dimension fl39[) and the initial condition 



15 



11 a s (Mw) . . 22 a s {M w ) 

z+(M w ) = l + — — — — , z-{M w ) = l-— — — — . (46) 

In this way we obtain (after transforming back to the basis 0\ 2 which is more commonly 
used) at fj, = m b ~ 5 GeV: 

d(m b ) = +1.0952 (LO) - 0.0250 (exact NLO) - 0.0017, 

C 2 (m b ) = -0.2250 (LO) + 0.0575 (exact NLO) + 0.0043. (47) 



These equations should be read in analogy to Eq. (|42j): the first two numbers are the exact 
LO and NLO results, the last number indicates the summed contribution beyond NLO from 
the anomalous dimensions in the NNA approximation^]. We conclude that this contribution 

3 'ln the NNA approximation (essentially corresponding to the large 0a limit of QCD) one should for con- 
sistency always use the LO running coupling a s (fi), neglecting higher order corrections to the QCD beta 
function. Both of the last numbers in Eq. ( |47| ) were nevertheless obtained by solving Eq. (^) with the NLO 
running coupling (and the anomalous dimensions exact up to NLO), since we wished to estimate the sum of 
all contributions beyond NLO in the NNA approximation. 
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is completely negligible as compared to the NLO correction. The reason for this is that the 
NLO contribution itself is dominated by the large matching correction ( filf ) so that the small 
modifications due to fl43| ) in the running become unimportant. 

(ii) B° — B° mixing 

The effective Hamiltonian for B° — B° - mixing is |23] 



He// = ^ M W KM 2 zf^)0 + . (48) 

It contains only the operator + for which the NNA approximation worked very well even 
at NLO. This leads to the conjecture that NNA estimates of higher orders might be to some 
extent reliable here. The initial condition for the Wilson coefficient z^_ B (fi) is a complicated 



function of the top mass which can be found in the appendix of |23|. Evaluating this function 
at m t (Mw) — 177 GeV, we obtain 

zf(M w ) = 2.624 + 2.798 a ^ Mw \ (49) 



The evolution down to the low scale with Eq. ( |4q ) then yields 

z BB {m b ) = 2.284 (LO) + 0.028 (exact NLO) + 0.006. (50) 

Again, the summed contributions beyond NLO are minute and can safely be neglected. 

From Eqs. ( f47| ) and ( |50| ) we conclude that the higher-order fermion-bubble contributions 
to the anomalous dimensions 7± are irrelevant in phenomenological applications. The reason 
for this is not the generic smallness of these additional contributions, cf. Eq. (|4"3|), but the 
smallness of the running in the current-current sector in general. Our analysis indicates that 
the termination of the perturbative series at NLO is in this case a very good approximation, 
with an associated theoretical error of almost zero. 

However, a word of caution is in order here. We examined in this paper only the large 
Nf contribution to the anomalous dimension. The result Eq. (|39| ) was a perfectly convergent 
series, essentially a product of gamma functions with a s Nf/6ir as argument, without any 
renormalon ambiguity. This is in fact the general structure of renormalization group functions 
in the large Nf limit 13, 19, 2^. Conversely, this implies that there is no n! growth of the 



perturbative coefficients and so there is no longer a good reason to believe that at sufficiently 
high orders the large Nf contributions become dominant. In the best case we can consider 
them as a guess of what may happen at higher orders. 

The n! growth is indeed not present (at least in dimensional regularization) in the divergent 
parts of Feynman diagrams but only in the finite parts - precisely the pieces we neglected in 
this calculation. This somewhat surprising statement can be verified directly from Eq. (^6|). 
The finite parts manifest themselves in corrections to the matching of the effective onto the full 
theory, i. e. they contribute to the initial conditions Ci(Mw) of the Wilson coefficients at the 
electroweak scale. So this initial condition is the place where one really expects renormalon 
singularities in the Borel plane. Hence, a complete large Nf analysis should certainly include 
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these matching corrections - a task which remains to be done in the future. If the general 
perturbative expansion of a Wilson coefficient is written as 

n=0 m=0 x * ' 

with expansion coefficients that are by themselves power series in Nf, 

C nm = J2cWN), (52) 

fc=0 

then the calculation of the anomalous dimension provides us, in principle, with all the coef- 
ficients C™" 1 with m > 0, i. e. with all the logarithmic pieces. The matching calculation on 
the other hand gives us the coefficients C^q 1 which will grow as (n — 1)! and will therefore 
presumably dominate the whole perturbative expansion at sufficiently high orders. Strictly 
speaking, this means that the whole renormalization group evolution from \i ~ Mw down to 
jU ~ mi,, i. e. the summation of large logarithms, becomes meaningless in the large Nf limit, 
since the perturbative series in this limit is (again: at sufficiently high orders) dominated by 
the constant and not by the logarithmic pieces, even if the latter are large. More work needs 
to be done in order to clarify these questions completely. 



6 Summary 

In this paper the anomalous dimensions of the current-current four-quark operators 0\ and 
2 were calculated in the large Nf limit in the MS (or MS) scheme. The calculation 
was done via direct computation of Feynman diagrams in dimensional regularization, with 
a special emphasis on the role of evanescent operators. Our findings agree with the a 2 s Nf 
part of the exact NLO result and may furthermore serve as an independent check for future 
NNLO calculations. Using the NNA approximation, we estimated with our result the total 
contribution of all orders beyond NLO to the anomalous dimension. This contribution was 
found to be of the same order of magnitude as the exact NLO correction itself. We also 
investigated to which extent this contribution affects physical quantities. As an example, we 
considered the Wilson coefficients relevant for nonleptonic B decays and B° — B° mixing. 
In these two cases, the effects were negligibly small and we concluded that the truncation of 
the perturbative series at NLO is a very good approximation in the current-current sector. 
Finally, future steps were discussed that will be necessary for a complete large Nf analysis of 
effective weak Hamiltonians. 
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Appendix: Evanescent operators 



As we observed in Sec. 3, evanescent structures as e. g. 

£i = 7*VV(1 - 75) ® 7m7^7a(1 - 75) - 4 (4 - e - e 2 ) -f(l - 75) ® 7 M (1 " 7s) (53) 

unavoidably arise while calculating matrix elements of weak four-quark operators^} For a 
consistent calculation these operators must be included from the very beginning,, i.e. one 
should start such a calculation from a "true" effctive Hamiltonian 

00 

H eff = C + J2Ce j £ j , (54) 

where the second term constitutes an infinite tower of evanescent operators vanishing in the 
limit e — > 0. To avoid an embarrassing amount of indices, here and in the following two 
simplifying assumptions are made: (i) we are dealing with only one physical operator 0, (ii) 
only one new evanescent operator arises at each loop, i. e. the n loop matrix elements can be 
written as 

n 

<0>W = a<0>W + 2><£,><°\ 

3=1 

00 

= Ci(C»)( )+^^(^) (0) , (55) 

where the sum in the first equation stops at n and not at infinity, and with some coefficients 
a, bi, Ci, dij which are computed by explicit calculation. Of course both assumptions can 
easily be relaxed. 

Now the crucial point is that we wish to avoid finite contributions from evanescent oper- 
ators to the physical matrix element 

00 

(Hf ff ) = C (0 R ) + J2Cs j (^)- (56) 
i=i 



In principle, this can be achieved by choosing the multiplicative operator renormalization 

(57) 




Zoo Zoe 
Zs.o Zeis 




in such a way that the vanish. Plugging Eq. (|57|) into Eq. (1561) and using the matrix ele- 



ments (j55|), this condition (together with the finiteness of {T~i^rf)) enables us to determine the 



L eff/ 

renormalization constant Z. Note that Z is the inverse of the usual operator renormalization 
constant Z. 



4 'lt is essential for the following argument to define the evanescent operators in such a way that Si _L O in 
D = 4 — 2e dimensions so that the coefficients Ci in Eq. (JHB]) start at 0(a s ) and not at 0(a®). This condition 
fixes the coefficient 4 (4 — e — e 2 ) in Eq. ( p3[ ) uniquely. In general, however, the O(e) part of this coefficient can 
be chosen arbitrarily (implying different definitions of the evanescent opertors) without affecting the physical 
results. For a detailed discussion on the subject of evanescent operators we refer to the work of Herrlich and 
Nierste @. 
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However, this procedure involves of course a finite renormalization in the Z^q in order 
to remove the finite contribution Ci{0)^ of the evanescents to the matrix element (55). The 
general structure of the renormalization constant Z at n loops is therefore 



V 







(58) 



where s means "only singular parts" and / "also finite parts". Note that finite O(e ) contri- 
butions to the bi and dij in Eq. (j55|) nevertheless yield no finite contribution to the matrix 
element due to the vanishing of the evanescents in four dimensions. This is the reason why 
the Zo£j an d ZEiSj m (§) contain only singular pieces. Furthermore, Zo£j = for j > n 
since only the first n evanescent operators contribute to the n loop matrix element (0)^ n \ 

However, there is a price to pay for the finite renormalization in (]58|). Namely, the finite 
parts of the Z^.q will in general result in an additional finite contribution to the anomalous 
dimension "foo of the physical operator O. Hence, the running and thereby the numerical 
value of the Wilson coefficient Cq in Eq. ( |56| ) is modified in a calculable way although the 
evanescent matrix elements themselves do not contribute. This subtle effect was first observed 
in this context in [jTj| . 

The relation between the renormalization constant Z = Z~ l and the corresponding 
anomalous dimension 
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~ dZ 

Zll — 

d[i 



(59) 



can easily be calculated at n loops. From Eq. (58) it follows for the general structure of Z 



( s 



\ 



\ 



(60) 



Using (p8|), (pOD and the finiteness of 7 one obtains in the usual way 



dZ, 



loo 



-9- 



(1) 

00 



dg 



n 
k=l 



7(1) ri dZ £ k O 
dg ■ 



(61) 



where we have expanded Z and Z in powers of 1/e, 
Z = z(°) + -Z( 1 )+0(e- 2 ), 



e 

1 Al 



Z = z^ + -Z^+0(e 



-2\ 



(62) 



The second term in Eq. (|6l]) would vanish without the finite contribution due to the evanes- 
cents, so Eq. (pl|) indeed makes sense as generalization of the usual formula (|T6|). 
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Now one can immediately conclude from Eq. ( |6l|) that there is no contribution to 700 
from the evanescents in the large Nf limit. This is for the following reasons: Nf enters the 
matrix elements ( |55| ) and therefore the renormahzation constants Z, Z only at 0(a 2 ), i.e. 

Ze k o = a s A k + a 2 s N f B k + 

Zo£ k = a s C k + a 2 s N f D k + .... (63) 

Multiplying these we see that the Nf terms in the evanescent contributions to the r.h.s. of 
Eq. ( |6l"| ) start with a^Nf, so they are suppressed compared to the a^Nj terms contained in 
Zqq. In the limit Nf — > 00 they can therefore be neglected. Note that this statement is only 
true if there are no 0(aP s ) terms present in Eq. (|63|) . In order to ensure this we defined the 
evanescent operators as described in the footnote above. 
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